MSUHEP-030303 



X 



Inclusive Jet Production, Parton Distributions, 
and the Search for New Physics 

Daniel Stump, Joey Huston, Jon Pumplin, and Wu-Ki Tung 

Department of Physics and Astronomy 
Michigan State University 
East Lansing, MI 48824 



cn 

o 
o 

' Department of Physics and Astronomy 



H.L. Lai 



Michigan State University 
East Lansing, MI 48824 and 



>• ' Department of Science Education 

m 

^ ^ Taipei Municipal Teachers College, Taiwan 

cn 
o 

cn 

o 

^7^. Argonne, IL 60439 

D 

^ ■ J. F. Owens 



Steve Kuhlmann 
Argonne National Laboratory 



Department of Physics 



^ \ Florida State University 

Tallahassee, FL 32306 
(Dated: March 3, 2003) 



1 



Abstract 

Jet production at the Tevatron probes some of the smallest distance scales currently accessible. 
A gluon distribution that is enhanced at large x compared to previous determinations provides 
a better description of the Run lb jet data from both CDF and D0. However, considerable 
uncertainty still remains regarding the gluon distribution at high x. In this paper, we examine the 
effects of this uncertainty, and of the remaining uncertainties in the NLO QCD theory, on jet cross 
section comparisons to Run lb data. We also calculate the range of contributions still possible 
from any new physics. Predictions are also made for the expanded kinematic range expected for 
the ongoing Run 2 at the Tevatron and for the LHC. 
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I. INTRODUCTION: JET CROSS SECTIONS AT THE TEVATRON 



Jet production at the Tevatron probes the highest momentum transfer region currently 
accessible. As this region is potentially sensitive to a wide variety of new physics, there was 
great interest when the inclusive jet cross section measured by the CDF collaboration in Run 
lb exhibited an excess in the highest Et range, when compared to NLO predictions using 



then-current parton distribution functions In an attempt to determine if conventional 
physics could explain the deviation, the CTEQ PDF fitting group carried out a global 
analysis using information from deep-inelastic scattering and Drell-Yan data, as well as 
from jet data from the Tevatron, but giving a higher statistical emphasis to the high Et jet 

n 

data from CDF [2j. NLO predictions using the resulting fit (CTEQ4HJ) reduced the size of 
the excess observed by CDF. The jet data from both CDF and D0 were also used in a more 
conventional fit (CTEQ4M) where no special emphasis was given to the high E^ data P|. 

The dominant subprocess that contributes to jet production at high Et is quark- 
(anti)quark scattering. However, the quark distributions in the corresponding x range are 
very well constrained by the precise DIS and DY data used in the global fits. Only the 
gluon distribution has the flexibility to change significantly in the high x region, and indeed 
the gluon distribution increases by roughly a factor of 2 for x values of about 0.5 in the 
CTEQ4HJ fit. The gluon-quark scattering subprocesses increase from approximately 20% 
of the total jet cross section at high Et using CTEQ4M to 40% using CTEQ4HJ 

The next group of PDF's from CTEQ (CTEQ5) l^] also contained two sets: CTEQ5M, 
the standard lowest solution, and CTEQ5HJ, defined by a similar statistical enhancement 
applied to the high Et data from CDF jsf and D0 . 

More recently, D0 has measured the inclusive jet cross section as a function of rapidity y 
over the range < |y| < 3 0]. This data set, comprised of 90 data points, has a greater sta- 
tistical power in the global fits than the CDF and D0 jet cross section measurements in the 
central region alone. The CTEQ6M fit j^l utilizes these D0 jet cross section measurements, 
along with the CDF measurements in the central rapidity region, as well as the most recent 
DIS data from HERA and existing fixed target DIS and DY data. An important step was 
taken to quantify the uncertainties of the parton distributions and their physical predictions 
with the inclusion of 40 eigenvector basis PDF sets, along with the best fit, to characterize 
the uncertainties in the neighborhood of the global minimum, using the Hessian method. 
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This paper is devoted to a detailed and more focused study of the uncertainties of the 
inclusive jet cross section due to global constraints, both at current energies, and more im- 
portantly for further QCD studies and the search for new physics at future collider programs. 
During the course of this work, we have produced an improved version of the CTEQ6 PDF's. 
The new minimum set, called CTEQ6.1M henceforth, provides a global fit that is almost 
equivalent in every respect to the published CTEQ6M ^, although some parton distribu- 
tions (e.g., the gluon) may deviate from CTEQ6M in some kinematic ranges by amounts 
that are well within the specified uncertainties. The more significant improvements are asso- 
ciated with some of the 40 eigenvector sets, which are made more symmetrical and reliable 
in CTEQ6.1M1 

The comparison between the D0 data for the 5 different rapidity intervals, and the predic- 
tions using the CTEQ6.1M PDF's, is shown in Fig.^ The theory predictions for CTEQ6M, 
CTEQ5M and CTEQ5HJ parton distributions are also shown. The greater statistical power 
that results from including the new D0 jet data, and the fact that all of the rapidity inter- 
vals prefer a larger gluon at high x, results in the increases in the cross sections found using 
NLO predictions based on CTEQ6M or CTEQ6.1M rather than CTEQ5M. 

This point deserves further emphasis. It is crucial to measure the jet cross section over a 
wide rapidity range: any new physics should contribute mostly to the central rapidity region, 
while a PDF explanation will contribute to all rapidity ranges. This will be addressed in 
more detail below. 

As discussed in Ref. Q], it is important to have a parameterization for the PDF's that 
is flexible enough to allow an increase in the high x gluon, without modifying too much 
the shape of the gluon distribution at lower values of x. A comparison of the CTEQ5M, 
CTEQ5HJ and CTEQ6M gluon distributions is shown in Fig.|^ As can be observed, the 
CTEQ6M gluon at high x lies between the gluons in CTEQ5M and in CTEQ5HJ, just as 
the CTEQ6M predictions for the inclusive jet cross sections at high Ej- lie between the 
predictions using the two other PDF's. The CTEQ6.1M gluon distribution is shown by 
dashed curves in Fig.|21 At the highest x values the function resembles CTEQ5HJ, and in 
the intermediate x range it resembles CTEQ6M. 



^ The CTEQ6.1M parton distributions will be available in the CTEQ table format at 
www.phys.psu.edu/ cteq/ and in the LHAPDF format at vircol.fnal.gov. 
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FIG. 1: The D0 inclusive jet cross section versus Et for five rapidity bins compared to NLO 
predictions using the CTEQ6.1M PDF's. The abscissa is Et in GeV. The ordinate is (data — 
theory) /theory. The curves show the CTEQ5M, CTEQ5HJ, and CTEQ6M predictions, as frac- 
tional differences compared to CTEQ6.1M. The error bars are the statistical and systematic errors 
combined in quadrature. 

The subprocess contributions to jet production are shown in the next figures, for the 
central rapidity bin (Fig.E)) and the forward rapidity bin (Fig.0)), for the parton distributions 
CTEQ5M and CTEQ6M. 

The dominant difference between the CTEQ5M and CTEQ6M sets is the increased gluon 
distribution in CTEQ6M at large values of x; the quark distributions are nearly unchanged. 
This shows up as a decrease in the fractional contributions of the quark-quark subprocesses 
and as increases in the fractional contributions of the quark-gluon and gluon-gluon subpro- 
cesses. It is interesting to note that the gq scattering subprocess is even more important at 
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FIG. 2: Left: The CTEQ5M, CTEQ5HJ, and CTEQ6M gluon distributions at = lOGeV^. 
Right: The ratios of CTEQ5M and CTEQ5HJ gluon distribution to that of CTEQ6M. The dashed 
curves show the CTEQ6.1M gluon distribution. 

high Et for the highest rapidity bin than for the central rapidity bin. 

The cross section for production of high-E'T- jets at the Tevatron depends on the parton 
distributions at large x. The relevant range of x may be estimated by considering the 
leading-order 2-^2 kinematics. Let xi and X2 be the momentum fractions of the two 
incoming partons. The outgoing parton 4-momenta are 

P'l = {pt cosh y I, pt, Pt sinh y i) , (1.1) 
pT = (pt cosh ^2, -PT,PTsinh 1/2) ; (1-2) 

the outgoing partons have rapidities yi and y2, and transverse momenta pt and —pr, respec- 
tively. If the first parton is identified with the observed jet then the jet rapidity is yj = yi 
and its transverse momentum is px- The 2 — > 2 kinematics determines the relation between 
the momentum fractions and outgoing parton variables, 

xi = ^ (e» + e^2) , (1.3) 
X2 = ^ ie-y^ + e-y^) . (1.4) 

These parametric equations yield a curve in the (xi, X2) plane, parametrized by the rapidity 
?/2 of the second parton, for the specified jet variables. 

Figure El shows the (xi,X2) values for central rapidity (?/j = 0) on the left and forward 
rapidity {ijj = 2) on the right. In the central case, the solid curve is the locus of {xi,X2) points 
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FIG. 3: The subprocess contributions to central jet production. 

for ■pT — 400 GeV and the dashed curve is for pr = 200 GeV. Both x\ and xi must be > 0.25 
to produce the high-pr jet. In the forward case the sohd curve corresponds to pt = 200 GeV 
and the dashed curve to Pt = 100 GeV. Here one parton must have large momentum fraction 
to produce the jet at large pr and forward rapidity. The second parton can have a low value 
of leading to the second jet being close in rapidity to the first jet, or a larger value of x, 
leading to configurations where the two jets are on opposite sides of the detector. Consider 
the contribution of the quark-gluon scattering subprocesses. For the configuration where 
one X is small and the other is large, the smaller x value most likely corresponds to a gluon 
and the larger to a quark. In this case the difference between the CTEQ5M and CTEQ6M 
distributions is small. On the other hand, for the configuration where both x values are 
large the increased gluon in CTEQ6M will enhance the cross section. The conclusion is 
that the increase in the cross section when going from CTEQ5M to CTEQ6M occurs in 
configurations where the forward jet is balanced in pt by a jet with roughly the opposite 
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FIG. 4: The subprocess contributions to jet production in the forward rapidity region, 
rapidity. 

Next, consider the question of where in the jet phase space one is most hkely to observe 
signs of compositeness or other new physics. The angular distributions of the dominant QCD 
subprocesses are sharply peaked at small scattering angles, whereas the models usually used 
to estimate compositeness have much flatter distributions. This means that the signal/QCD 
ratio will be largest for 90° scattering in the parton-parton rest frame. In the overall hadron- 
hadron frame, this corresponds both jets having the same rapidity. Using the lowest order 
kinematics described previously, the squared dijet mass is given by 

Mj = 2p|[l + cosh(yi-y2)]. 

Setting yi — 1/2 yields Mjj — 4p^. This shows that maximizing the dijet mass means maxi- 
mizing pt- Since the pr reach is maximized in the central region, this is the region where the 
highest dijet masses can be reached with the lowest background for new physics signals (at 
least those with a flat angular distribution) . The model discussed for compositeness later in 
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this paper gives an explicit example of this behavior. On the other hand, as shown by the 
kinematic examples in Figs. El modifications of the gluon distribution at high values of x 
will affect the high-pj^ jet cross section at all rapidities. 
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FIG. 5: Left: Parton momentum fractions xi and X2 for yj = 0; the solid and dashed curves are 
for Et = 400 and 200 GeV, respectively. Right: Parton momentum fractions xi and X2 for yj = 2; 
the solid and dashed curves are for Et = 200 and 100 GeV, respectively. 
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II. UNCERTAINTY ANALYSIS OF THE CROSS SECTION FOR INCLUSIVE 
JET PRODUCTION AT THE TEVATRON 



As discussed in the previous section, the cross section for inclusive jet production has 
been measured by the CDF and D0 collaborations at the Tevatron. These results have 
been used in the latest CTEQ global analysis of PDF's, leading to the CTEQ6 (and now 
6.1M) parton distributions. In this section we analyze the uncertainty of the theoretical 
QCD cross section due to PDF uncertainty. 



A. The Hessian method of uncertainty analysis 

The Hessian method of uncertainty analysis has been described in detail previously [3| ■ 
For completeness, we summarize the method briefly. The parton distributions are con- 
structed by a method of chi-square minimization with fitting of systematic errors. A chi- 
square function x'^ is defined by 

X'' = Y.^lia,r) (2.1) 

e 

,) = J2 "'■'^^''f' + E rl (2.2) 

i » k 

where e labels an experimental data set and i labels a data point in the data set. Di 
is the data value, is the uncorrelated error, and jSki is the kth correlated systematic 
error; these numbers are published by the experimental collaboration. Then Tj(a) is the 
theoretical value, a function of a set of n PDF parameters, {ai, . . . , a„}. Also, {r^} is a set 
of Gaussian random variables and r^Pki is a (correlated) shift applied to to represent 
the kth systematic error. We minimize the function x'^('j^) ^) with respect to both the PDF 
parameters {a} and the systematic shift variables {r^}. The result yields both the standard 
PDF model with parameters {oq}, and the optimal shifts {ffe} to brine theory and data into 
agreement. This minimum of x'^ represents the best fit to the data 
The Hessian matrix is 

1 92 £2 

Hij = 2 gada- = ^(o))- (2-3) 

This matrix determines the behavior of x^{a) in the neighborhood of the minimum. We 
have developed an iterative method for calculating the Hessian accurately, leading to an 
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improvement in the minimization 12] . The point {ao} in the n-dimensional parameter 
space — where x^{a) is minimum — is the best fit to the global data set. However, points in 
some small neighborhood of {ao} are also acceptable fits. The variation within the acceptable 
neighborhood, of a PDF or a physical prediction, is our measure of the uncertainty of the 
quantity. 

In the Hessian method of uncertainty analysis, we first define n special directions in the 
parameter space, namely the n eigenvectors of the Hessian. Then for each eigenvector we 
have two displacements from {oq} (in the + and — directions along the vector) denoted 
{af} and {a^} for the ith eigenvector. At these points, = Xo + where Xo = X^{(^o) 
= the minimum, and T is a parameter called the tolerance. We consider any PDF set 
with — Xo < T'^ to be an acceptable fit to the global data set. In particular, the 2n 
PDF sets {af}, which are called the eigenvector basis sets, span the parameter space in the 
neighborhood of the minimum. 

The appropriate choice of tolerance T cannot be decided without a further, more detailed, 
analysis of the quality of the global fits. After studying a number of examples [3, ll^ we 
conclude that a rather large tolerance, T ~ 10, represents a realistic estimate of the PDF 
uncertainty.^ 

Any quantity X that depends on PDF's has a predicted value Xq = X(ao) and a range 
of uncertainty SX. A simple measure of 5X is the spread of the values of X{af) for the 2n 
eigenvector basis sets. However, a more complete uncertainty range is between the minimum 
and maximum values of X for all points with x^ — Xo < T"^. It can be shown that in a linear 
approximation these bounds are Xq ± 6X where 

(.X)=^T^^(i/-),|f^; (2.4) 

or, in terms of the eigenvector basis sets, 

{SXf = -Y,[X{ar)-X{ar)]\ (2.5) 
fc=i 

We refer to ()2.5|1 as the Master Equation for the uncertainty of X in the Hessian method [3| • 
Equation (j2.5p is based on a linear approximation: x^{a) is assumed to be a quadratic 



^ Different choices of the tolerance have been made by other groups, e.g., T ~ y/AO by MRST [lo| . 
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function of the parameters {a}, and X{a) is assumed to be linear. We find that these 
approximations are not strictly valid, so instead we calculate asymmetric bounds by 

n 

5X^ = J2[X{at)-X{ao)Y; (2.6) 

k=l 

the range of uncertainty of X is {Xq — 6X^,Xo + (In ()2.6|) {af} and {a^} are the 

displaced points where X > Xq and X < Xq, respectively.) 

The uncertainties on the gluon and u-quark distributions are shown in Fig.lHl for = 
10 GeV^. The w-quark distribution is tightly constrained for x < 0.8, whereas the uncertainty 
on the gluon distribution grows to a factor of 2 or larger for x values greater than ~ 0.4. 
For comparison, the CTEQ5M, CTEQ5HJ, MRST2001 Q and MRST2002 10] PDF's are 
shown, compared to CTEQ6.1M, in Fig.H for = lOGeV^. 



Gluon uncertainty U-quark uncertainty 
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FIG. 6: PDF uncertainties. The curves are the ratios of the 40 eigenvector basis sets to the standard 
set (CTEQ6.1M). The largest variations of the gluon function are the + and — displacements along 
eigenvector 15. 



B. The Tevatron jet cross section 

To illustrate our method, and the importance of the systematic errors, Fig.|Hl shows 
the fractional difference between the data obtained by CDF in Run lb, and the theory 
with CTEQ6.1M partons. The error bars are the statistical errors only. The left panel is 
{Di — Ti)/Ti. Note that there are some systematic differences between the CDF jet data and 
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FIG. 7: Comparison of several standard PDF sets. 



the theory, both in terms of normahzation and in shape. ^ The right panel shows {Di — Ti)/Ti 
where Di is the shifted data, 

n 

D, = D,-J2rAf, (2.7) 

k=l 

i.e., the optimal systematic shifts determined by the global fit have been subtracted from 
the data values. There is no systematic difference between Di and T,, as expected. The sizes 
of the systematic shifts are, as expected, comparable to the standard deviations published 
by the CDF collaboration; that is, is of order 1 for each of the systematic errors. 
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^ It is important to note all Run 1 CDF cross sections are intrinsically larger than comparable D0 cross 
sections by approximately 6% due to different assumptions regarding the total inelastic cross section at 
1.8 TeV lig. 
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FIG. 8: The left panel shows {Di — Ti)/Ti, the fractional difference between data and theory. The 
right panel shows {Di — Ti)/Ti where Di is the data with systematic shifts subtracted. The error 
bars are the statistical errors only. 

Figure ini shows the cross section for inclusive jet production at the Tevatron, integrated 
over the rapidity range 0.1 < \ri\ < 0.7, as a function of jet transverse energy Et-, for 
CTEQ6.1M partons and for the 40 eigenvector sets. (The rapidity interval is that for which 
the cross section has been published by the CDF collaboration.) The spread of the 41 curves 
in Fig.lHlis a simple estimate of the PDF uncertainty. The CDF data points are superimposed 
on the calculation. 
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FIG. 9: Calculations of da/dET for the 40 eigenvector sets. 

Figure UHl shows the fractional differences [T{af) — T(ao)] /T{ao) between the 40 sets 
and the standard prediction (CTEQ6.1M). The points are the CDF data compared to 
CTEQ6.1M; the error bars are statistical errors and the data points are plotted without 
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the systematic shifts. Also shown is the prediction using CTEQ5HJ partons. The range of 
uncertainty encompasses both the CDF data and the CTEQ5HJ predictions. 
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FIG. 10: Calculations for the 40 eigenvector basis sets, plotted as fractional differences compared 
to CTEQ6.1M. The points are the CDF measurements. The error bars are statistical errors only; 
the systematic shifts are not subtracted from the data. The dashed curve is CTEQ5HJ compared 
to CTEQ6.1M. 

Figure ^2 shows the 40 eigenvector basis sets separately. For each eigenvector the two 
curves are the positive and negative displacements. The Hessian method assumes that 
the variations are approximately linear. If the linear approximation were strictly valid the 
curves would be mirror images about zero. Fi gure ITT] implies that the linear approximation is 
reasonably accurate for most eigenvectors, especially those corresponding to the ten largest 
eigenvalues (1 - 10). 

Figure ^1 shows the overall PDF uncertainty calculated from the "Master Equation," 
asymmetrically for positive and negative differences. The ordinate is the fractional difference 
between the extreme cross section and the central (CTEQ6.1M) prediction. The outer curves 
are the upper and lower bounds of the uncertainty band. The CDF data is superimposed. 

The overall uncertainty in the calculation of {da /dEx) cbf is dominated by one of the 
eigenvector directions. The dominant eigenvector is "eigenvector 15," i.e., that for which 
the eigenvalue of the Hessian is 15th in order of magnitude. The + and — displacements in 
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FIG. 11: The cross sections for + and — displacements along the 20 eigenvectors, plotted as 
fractional differences compared to CTEQ6.1M. 

parameter space along this direction have a much larger effect on the value of the cross section 
than any of the other eigenvector directions. In other words, {da / dEx) cdf is most sensitive 
to the parameter variations along this eigenvector. Figure El shows the gluon distribution 
for CTEQ6.1M and for the eigenvector basis sets in the + and — directions along eigenvector 
15. Eigenvector 15, like the 19 others, corresponds to some linear combination of all of the 
20 fitting parameters. The exceptional feature of eigenvector 15 is that it represents the 
largest change in the high-x gluon behavior. 

In the new CTEQ6.1M analysis, the eigenvector basis sets have been carefully constructed 
to have = Xo + 100.^ But is the global chi-square function. The criterion Ax^ = 100 is 
only a reasonable rule of thumb, and a more thorough analysis — inspecting the fit for each 
data set — is needed to justify the extreme sets a^^ and a^^. We must ask whether these 

^ The earlier CTEQ6 analysis made an approximation of the NLO cross sections — holding the X-factor 
tables for DIS and DY processes constant as the PDF parameters varied in the neighborhood of the 
minimum — such that Ax^ was not strictly constant at for the Hessian PDF sets. The CTEQ6.1M 
analysis is preferred although the changes from CTEQ6 are small. 
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FIG. 12: The overall PDF uncertainty of the inclusive jet cross section is plotted as the two red 
curves. The points are CDF data; the error bars are statistical only and the systematic shifts have 
not been subtracted. 

variations could be ruled out by any specific experiment. 

The theoretical cross section most affected by displacements along eigenvector 15 is the 
inclusive jet cross section. Table IJ shows the values for the CDF and D0 data for three 
PDF sets ( ao, and 015.) Because of the significant systematic uncertainties in the exper- 
imental data, represented by the error correlation matrices published by the collaborations, 
we might conclude that the extreme sets af^ and a^^ are still acceptable fits to the data 
within the systematic errors. 



PDF set 


D0 


D0 x^/N 


CDF 


CDF xV^ 






{N = 90) 




{N = 33) 


CTEQ6.1M 


52.26 


0.581 


57.67 


1.748 


EV set 15+ 


111.90 


1.243 


53.05 


1.608 


EV set 15- 


39.17 


0.435 


67.35 


2.041 



TABLE I: x^ /N for the standard PDF's (CTEQ6.1M) and the extreme variations in the + and — 
directions along eigenvector 15. 

As implied by the table above, much of the increase in (50 out of 100) for eigenvector 
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FIG. 13: Eigenvector 15. Comparison between the gluon distributions for the + and — displace- 
ments along eignevector 15 and the standard set CTEQ6.1M. 

+15 comes from the D0 jet cross section. The Run lb D0 cross section can tolerate 
some increase compared to the CTEQ6.1M predictions, but the increase corresponding to 
eigenvector +15 is in fact too extreme. It is not surprising that such a large portion of the 
increase in comes from the D0 jet data: as shown in Fig.^J this eigenvector probes 
primarily the high-x gluon distribution and the other data sets (with the exception of the 
CDF jet data) are fairly insensitive to this direction. The opposite direction, eigenvector 
— 15, gives the lowest fo^^ the D0 jet data. (Of course the x^'s of other experiments are 
larger so that the net change is 100.) The cross section in this case is less than that for 
CTEQ6.1M, by as much as 10 - 20%. The fact that x^/N is much less than 1 for the D0 
data suggests that the size and correlation of systematic errors are not accurately understood 
and thus CTEQ6.1M gives a satisfactory fit although not the best fit. 
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Figures for the D0 jet data analogous to Figs. M and HUl are shown in Figs. [T31 and IT31 
Figure IT^ shows the inclusive jet cross section as a function of Et for the five rapidity bins of 
the D0 measurements. Calculations for all 40 eigenvector basis sets are superimposed. The 
points are the D0 data and the error bars equal the single-point statistical and systematic 
errors added in quadrature.^ Figure ^1 shows the inclusive jet cross section as a function 
of Et for the five rapidity bins, plotted as fractional differences compared to the central 
(CTEQ6.1M) model. 
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Et [GeV] 



^ The systematic errors for the D0 data are supphed by the collaboration only in the form of correlation 
matrices. Unlike the CDF data, for which the standard deviations of individual systematic errors have 
been published, we are unable to separate the optimal systematic shifts for subtraction from the D0 data. 
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FIG. 14: The inclusive jet cross section as a function of Et for the five rapidity bins of the Run lb 
D0 measurements. Predictions of all 40 eigenvector basis sets are superimposed. The points are 
the D0 data, and the error bars are the statistical and systematic errors combined in quadrature. 
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FIG. 15: The D0 inclusive jet cross section as a function of Et for five rapidity bins, plotted as 
fractional differences compared to the central (CTEQ6.1M) model. The abscissa is Et in GeV and 
the ordinate is [T{a) — T{ao)] /T(ao) for the curves and {D — T)/T for the points. 



21 



III. RELIABILITY OF NEXT-TO-LEADING-ORDER PREDICTIONS 



Thus far, the discussion of uncertainties has dealt primarily with the issue of the prop- 
agation of experimental errors on the observables and their effect on the parameters of the 
parton distributions. There are, of course, issues related to theoretical uncertainties both in 
the global fit itself and in the theoretical calculations for jet production. However, it is far 
more difficult to quantify these. We have previously described studies of some theoretical 
uncertainties |8i]. For example, we have chosen cuts on Q and W for the DIS data sets 
used for CTEQ6 and earlier analyses in order to reduce the theoretical uncertainties due to 
effects from higher twist and related terms. Also, we have varied the functional form of the 
parametrization to ensure that sufficient flexibility to cover the wide kinematic range of x 
and Q is available. In this section we concentrate on the theoretical uncertainties related 
to NLO calculations for jet production. The response of the predictions to variations in the 
renormalization and factorizations scales yields some information on the theoretical uncer- 
tainties. In this section the scale dependence of the jet cross sections and the effects on the 
fits of varying the scales will be discussed. Some estimates of the possible effects of two-loop 
corrections will also be presented. 

A. Scale dependence of the jet cross section 

Leading-logarithm calculations of the jet cross section at high-ii^j' generally decrease 
monotonically as the factorization and renormalization scales are increased. The scale vari- 
ation of next-to-leading-logarithm calculations is generally reduced as compared to that of 
the leading-log calculations. A brief review of how this comes about is presented in the 
Appendix. 

The increased rapidity coverage available with the recent D0 data has provided addi- 
tional constraints on the gluon distribution. However, if the theoretical uncertainties become 
large at high rapidity, then these constraints may not be as useful as first believed. As one 
approaches the edge of phase space, it might be anticipated that there will be large loga- 
rithms due to the constraints on multiple gluon emission. The effects of these large logs will 
be investigated in Sec. (jlll C|) . Their effects may show up already in the next-to-leading-log 
calculations used in the global fits, so it is important to look for any anomalous behavior in 
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FIG. 16: Comparison of the scale dependence of the jet cross section at Et=70 GeV/c for the 
leading- and next-to-leading-logarithm calculations 

the predictions. 

The first step is to examine the scale dependence of the jet cross section predictions in 
the highest rapidity bin of the D0 data. For ease of presentation, the factorization and 
renormahzation scales have been chosen to be the same. In the jet cross section calculations 
used for the CTEQ PDF fits, and in this paper, a scale = Hf = Et/I has been used. 
Fig. [TBI shows the comparison of the scale dependences of the leading-log and next-to-leading- 
log calculations at Et = 70GeV/c. The monotonically decreasing behavior of the leading- 
log results is evident, as is the roughly quadratic dependence of the next-to-leading-log 
curve. Both behaviors are consistent with the expectations discussed in the Appendix. 
Similar results are shown in Fig.^lfor E^ = 170GeV/c. Both figures show that the scale 
dependence of the next-to-leading-log calculation is greatly reduced with respect to that of 
the leading- log calculation, provided that one avoids scale choices much smaller that /i = Ex. 
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FIG. 17: Comparison of the scale dependence of the jet cross section at £'^=170 GeV/c for the 
leading- and next-to-leading-logarithm calculations 

It is interesting to note that the ii"-factor (NLO/LO) for fi = Et/2 is approximately 
unity for Et = 70GeV/c, while it is significantly less for Et = 170GeV/c. The reasons for 
this are discussed in more detail in Appendix A. 



B. Fits with different scales 



The next point to be addressed is the effect on the fits of varying the choice of scale. 
To investigate this point, the CTEQ6M fits^ were repeated with scale choices of /i = E^ 
and 2Et- As compared to the original fit with n = Et/2 the latter two yielded increases 

^ This procedure has not been repeated with the CTEQ6.1M PDF's. However, since the changes in the 
region of the D0 and CDF data were minimal, the scale dependence in the 6M and 6.1M sets should be 
essentially the same. 
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in chi-square of 16 and 70, respectively. Thus, the choice of Et/I used in the CTEQ6M 
fits yields the best chi-square of the three choices, but this does not constitute a fit for the 
optimal scale. Note, however, that the rapid increase in the scale dependence for choices 
below about Et/I shown in Figs.lTHl and 1171 precludes fits with significantly smaller scales 
while the larger scale choices show an increase in chi-square. 

The gluon distributions from these fits are shown in Fig.^J relative to the standard 
CTEQ6M distribution corresponding to fi = Et/2. The effects at values of x below about 
0.4 are minimal, but the effects at larger x values are significant. In particular, the reduction 
in the cross section due to the larger scale choice of 2Et is compensated by an increase of the 
gluon distribution. However, it should be noted that the chi-square increased by 70 units 
for this fit. This is still less than the chi-square increase of 100 that determines the 40 error 
PDF's. 

A comparison of Fig. ^1 to Fig. IHl indicates that the uncertainty on the gluon distribution 
due to the scale choice for jet production in the global fits is everywhere less than the 
uncertainty from the treatment of the experimental errors. 
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FIG. 18: Comparison of fitted gluon distributions for fits made with different scale choices for 
the jet cross section. Comparison to Fig.lHl shows that this theoretical error is smaher than the 
uncertainty due to experimental errors. 

It is important to bear in mind that the gluon distribution itself is not an observable 
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quantity. Rather, the gluon and quark distributions contribute via convolutions with ap- 
propriate hard scattering cross sections. Accordingly, the impression provided by Fig. 
may be somewhat misleading as much of the gluon variation is offset by changes in the 
hard scattering cross section with which it is convoluted, with the end result being that the 
variation in the physical observables is much less. 

The jet cross section predictions using scales of -Et/2, Et and 2Et are shown in Figs.lT^ 
and|2nifor the D0 rapidity regions of 0.0 - 0.5 and 2.0 - 3.0. These curves may be compared 
to those for bins 1 and 5 in Fig. ^| The scale dependence variations are generally within 
the bands formed by the 40 eigenvector basis sets, indicating that the uncertainty introduced 
by the scale dependence of the jet cross sections is less than or on the order of that coming 
from the experimental errors. 
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FIG. 19: Scale dependence of the jet cross section prediction for < y < 0.5. 



C. Threshold resummation 

As noted in the introduction to this section, there are potentially significant contributions 
to the jet cross resulting from phase space limitations on gluon emission as one goes to 
large and large rapidity. Such corrections can be treated by threshold resummation 
techniques In Ref. |l7| these corrections for jet production were estimated at the two- 
loop level for central rapidity values and found to be small. The calculation has been 
repeated for the largest rapidity bin covered in the D0 data and the results are shown 
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FIG. 20: Scale dependence of the jet cross section prediction for 2 < y < 3. 

in Figs.l^ and EH by the dashed Unes. As was the case in the central region, the two-loop 
estimates of the next-to-next-to-leading-logarithm contributions are not large for the scale 
choices considered in this analysis. 

To summarize the results of this section, the next-to-leading-logarithm jet calculations 
show the expected reduction in scale dependence, even for large values of Et in the highest 
rapidity bin in the existing data set. Nevertheless, the remaining scale variation does show 
up in an increased uncertainty in the gluon distribution at large values of x. This uncertainty 
is, of course, partially offset by the compensating scale dependence of the hard scattering 
cross section. Even so, the scale dependence of the next-to-leading-log calculations remains 
as a contributing factor for the overall uncertainty of the predictions for new processes. 
Finally, the next-to-next-to-leading-log corrections estimated by expanding the threshold 
resummation predictions for jet production remain small, even at high Et in the highest 
rapidity bin considered, lending support to the use of the next-to-leading-log expressions in 
the global fit. 
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FIG. 21: Same as Fig. with the addition of the estimated two-loop contribution expected from 
threshold resummation. 
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FIG. 22: Same as Fig. (|17j) with the addition of the estimated two-loop contribution expected from 
threshold resummation. 
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IV. ROOM FOR NEW PHYSICS IN RUN IB 



Hadron collider jet data have traditionally been used to constrain models of new physics, 
usually models of quark substructure. The measurements include dijet angular distributions 
that are insensitive to theoretical uncertainties due to parton distributions Inclusive jet 
cross sections or dijet mass distributions typically provide better probes of quark composite- 
ness, but are sensitive to parton distributions at large x which have large uncertainties 
Currently the best limits are provided by a hybrid approach that uses the ratio of dijets in 
different rapidity bins j^]. This approach has the advantage that many theoretical and ex- 
perimental systematic errors cancel out in the ratio. It has the disadvantage that it removes 
some useful information, such as the absolute values ot the cross sections in the two rapidity 
intervals. 

There are several reasons to revisit this issue. First, the published limits are based on 
parton distributions several generations old, and the best-fit gluon distribution at large 
X is significantly changed. Second, the recent jet measurements at large rapidity provide 
constraints on the large x partons while being largely insensitive to new physics. Finally, 
and most important, the global analysis machinery has never been used in these analyses to 
refit the parton distributions in the presence of a new component to the QCD calculation. If 
the constraints on the parton distributions are not sufficient, it is possible that new physics 
could hide inside the PDF uncertainties, or current limits may be weakened. 

Compositeness signals are typically parameterized by a mass scale A which characterizes 
the quark substructure coupling. The coupling is approximated by a four-Fermi contact 
interaction giving rise to an effective Lagrangian. On 
coupling of quarks and anti-quarks has been calculated 



H';he term describing left-handed 



Cg, = A{27T/A'){qLYqL){qLl,qL) where A = ±l. (4.1) 

The simplest squared amplitude is the quark anti-quark s-channel annihilation process, 
with the first term the normal QCD interaction and the second due to substructure: 

I Ainu dd) r= ^aim^^^ + §. (4.2) 

Since u is negative, constructive (destructive) interference in this process occurs when 
A = —1 (+1). This is generally true for the other scattering terms as well, which will 
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be included in the calculations. As in the experimental papers, we will use the quark com- 



positeness implementation in PYTHIA 



22| to form the ratio of QCD+substructure/QCD, 



then multiply this ratio by the NLO QCD calculation for each iteration of parton distribu- 
tion as the PDF parameters are being fit. As expected, we find larger cross sections with the 
constructive interference choice A = —1, approximately a factor of 10 larger. However, the 
ratio of the first two eta bins used by D0 is slightly larger with the destructive interference 
term A = 1, and this is what is used for the current best limit of A > 2.7 TeV [2^. Therefore 
we have used the destructive interference choice for illustration in the rest of this paper. 

Figure ESI shows the cross section for jet production including a contact interaction, for 
three values of A: 1.6, 2.0 and 2.4 TeV. The cross section is plotted as a ratio to the pure QCD 
prediction with the CTEQ6.1M parton distributions. The range of kinematic parameters 
are the same as the D0 data, and the D0 data is superimposed on the plots with error bars 
of statistical and systematic errors combined in quadrature. 

If the effect of the contact interaction is smaller than the PDF uncertainty of the QCD 
prediction then the data cannot rule out the compositeness model; any difference between 
theory and data could be explained either as a PDF effect or as a sign of new physics. 
However, if the compositeness model disagrees with the data by an amount that is larger 
than the PDF uncertainty, the model is ruled out. The PDF uncertainty of the QCD 
prediction is comparable to the experimental error bars (see Fig. [T3|) . so we may use those 
error bars to judge whether the effect of the contact interaction is smaller than the PDF 
uncertainty. The contact interaction is seen to be ruled out for A = 1.6 TeV, primarily from 
the disagreement with the central rapidity interval, < rj < 0.5. The cross section at large 
rapidity (77 > 1.5) is not sensitive to the contact interaction. 

The largest two A values, A = 2.0, 2.4 TeV, are not inconsistent with the D0 data on 
inclusive jet production because the effect of the contact interaction is within the PDF 
uncertainty. The global fit for A = 2.4 TeV is even slightly better than for the pure QCD 
model, although the difference in the overall quality of the fit {Axgiohni = 9) much less 
than our standard tolerance for PDF uncertainty (A^gj^j^j^j < 100). 

If there is a contact interaction, then the CTEQ6.1M PDF's would be inaccurate because 
the data were fit to a pure QCD model. Therefore a complete analysis of the contact interac- 
tion requires that the PDF's be refit, comparing the data to the theory with compositeness. 
Figure |S1 shows the final result of this analysis. The cross section for jet production with 
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FIG. 23: Comparison of three models in which a hypothetical contact interaction contributes to jet 
production. The ordinate is the fractional difference between the model and the standard theory; 
the abscissa is the jet Ej- in GeV. The contact interaction has A = 1.6, 2.0, and 2.4 TeV for the 
three curves. In all cases the PDF's are the standard CTEQ6.1M. The D0 data are superimposed 
for comparison. 

the D0 kinematic parameters is again shown for three values of the A parameter (1.6, 2.0 
and 2.4 TeV) as a ratio to the pure QCD model with CTEQ6.1M. In Fig.Elthe PDF's for 
each A value come from a global fit that includes the contact interaction in the theory of 
the jet cross section. The smallest value of A is now more consistent with the data, but still 
ruled out. 

The limit derived from this analysis, A > 1.6 TeV, is not as stron g as those from either 
the CDF or D0 analyses that use angular distribution information [l8|. This is a result 
of the remaining uncertainty in the PDFs, and the fact that in this model of new physics 
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A [TeV] 


PDF set 


global 


D0 x^/N 


CDF xViV 


1.6 


CTEQ6.1M 


2111 


0.939 


5.45 


2.0 


CTEQ6.1M 


1966 


0.566 


2.07 


2.4 


CTEQ6.1M 


1948 


0.493 


1.73 


oo 


CTEQ6.1M 


1957 


0.590 


1.68 


1.6 


best refit 


2079 


1.230 


3.63 


2.0 


best refit 


1965 


0.576 


2.00 


2.4 


best refit 


1947 


0.493 


1.73 



TABLE II: x^'s for three theoretical models with a contact interaction. 

the angular distributions are quite different than QCD. Otlier deviations from the Standard 
Model might change the absolute cross sections while maintaining a more QCD-like angular 
distribution. In this case this analysis is the most in-depth study so far of allowed deviations 
in jet cross sections due to new physics. 

The limit derived from this analysis is also less than that obtained by the D0 collaboration 
when using the ratio of the dijet cross sections in the two rapidity bins (0.0—0.5 and 0.5 — 1.0). 
However, it should be noted that at least part of the higher limit from D0 may be due to 
the dijet cross section for 0.5 < r/ < 1.0 being somewhat larger than the standard model 
prediction, while the dijet cross section for 0.0 < t] < 0.5 is in better agreement with the SM 
prediction. The larger the ratio of the 0.5 — 1.0 to 0.0 — 0.5 rapidity intervals, the stronger 
the compositeness limit. The detailed information about the absolute values of the cross 
sections in each rapidity bin has not been used. 

Table HTl shows how the x^'s are affected by the contact interaction. 
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FIG. 24: Comparison of three models in which a hypothetical contact interaction contributes to 
jet production, with refitting of the PDF's. The curves have the same meaning as in Fig. 1231 The 
PDF's for the three models have been refit to the global data set to give the best fit including the 
contact interaction in jet production. 
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V. PREDICTIONS FOR RUN 2 



The increase in the center-of-mass energy (from 1.8 to 1.96 TeV) and the increased 
luminosity expected in Run 2 (from approximately 100 pb^^ to 2-15fb^^ will result in a 
dramatically larger kinematic range for measuring jet production, and thus for searching for 
possible new physics as well. Also, in Run 2, the new CDF endplug calorimeter will allow 
for an extension of measurements of the inclusive jet cross section to the forward rapidity 
region, similar to the analysis that has been performed by D0 in Run lb. 
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FIG. 25: Run 2 cross section for the D0 rapidity bins. Left: The 40 eigenvector basis sets. Right: 
The central prediction and overall uncertainty band. 

The predictions for the Run 2 jet cross section as a function of transverse momentum 
(using the D0 rapidity intervals) are shown in Fig.|2Hl With 2fb~^ (Run 2a) a measurement 
of the inclusive jet cross section out to a transverse energy of 600 GeV is possible, while with 
15fb~^ (Run 2b) the measurement can be extended to almost 650 GeV. The predictions are 
displayed on the left graph in the form of the curves for the 40 eigenvector basis sets of 
CTEQ6.1M. The central prediction and overall uncertainty band from the Master Equation 
fl2.6|) are shown separately on the right graph. Figure EHl shows the uncertainty range 
of the Run 2 cross section by displaying the ratios of the eigenvector basis sets to the 
central prediction. Near the kinematic limit, the PDF uncertainties allow the cross section 
predictions to be as much as a factor of 2 larger than those from the central fit (CTEQ6.1M), 
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again due to eigenvector 15. 
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FIG. 26: Uncertainty range of the Run 2 cross section for the D0 rapidity bins. The curves show 
the ratios of the 40 eigenvector basis sets to the central (CTEQ6.1M) prediction (ordinate) versus 
Pt in GeV (abscissa). 



For the Run 2 jet cross section predictions, we have used the Run 2 cone jet algo- 



rithm 
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26| which uses 4-vector kinematics both to define the jet and to specify its trans- 
verse momentum and rapidity. This algorithm has been adopted by both the CDF and D0 
collaborations for their inclusive Run 2 jet results. Note that this algorithm specifies the 
transverse momentum pt rather than the transverse energy E'j^. 

For completeness we also show the plots of the predictions of the inclusive jet cross section 
at the Tevatron Run 2 for the CDF choice of rapidity bins. Figure EZI shows the cross section 



^ The Run 2 jet algorithm results in cross section predictions for CDF that are ~ 7% lower than those using 
the Run 1 algorithm, over most of the kinematic range. 
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FIG. 27: Run 2 cross section for the CDF rapidity bins. 

While the statistical errors in Runs 2a and 2b will be small, so that the jet cross section 
can be extended to high p^, the experience from Run lb shows that the systematic errors 
will have a dominant effect on the significance of the results. However, we can anticipate 
that the PDF's, especially the gluon distribution, will be pinned down at large x by these 
data. The experimental collaborations will make the data most useful for global analysis by 
providing detailed information on systematic uncertainties. 
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FIG. 28: Uncertainty range of the Run 2 cross section for the CDF rapidity bins. The curves 
show the ratios of the 40 eigenvector basis sets compared to the central (CTEQ6.1M) prediction 
(ordinate) versus pT in GeV (abscissa) . 

A. The ratio of Run 2 to Run lb cross sections 

As mentioned previously, the Run 2 center-of-mass energy ^/s is higher than that from 
Run lb. The increase in ^/s will have only a modest impact on most physics cross sections. 
However, it will lead to a rather large increase in the production of high-ii^T jets. It may be 
useful to examine the ratio of jet production (for the same Et values) at the two different 
energies. In such a ratio, many of the theoretical (and experimental) errors will cancel. 
The Run 2 to Run lb ratio for the 40 sets of PDF's is shown in Fig. EHl for the central 
rapidity range of CDF (0.1 < < 0.7). As can be observed the ratio has a rather narrow 
theoretical error band from PDF uncertainties, i.e., smaller uncertainty than that of the 
absolute prediction of the inclusive cross section.^ 



^ Here we have multiplied our Run 2 predictions by a factor of 1.07 so that the comparison effectively uses 
the Run 1 algorithm for both the numerator and denominator. This is done for convenience since CDF 
will present its first results for the Run 2 jet cross section using the Run 1 algorithm. 
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FIG. 29: The ratio of the Run 2 to Run lb cross sections for the central rapidity range (0.1 < 
\ri\ < 0.7). The Run 2 cross section has been multiplied by 1.07 to take into account the difference 
in the jet definitions for the Run 2 and Run lb analyses. (The small fluctuations are due to the 
Monte Carlo errors in the calculations of the K-factors (NLO/LO) at each energy.) 
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VI. INCLUSIVE JET PRODUCTION AT THE LHC 



The increase of the center-of-mass energy to 14 TeV at the LHC will result in a dramat- 
ically larger accessible kinematic range. Inclusive jet cross sections will be measured out to 
transverse momentum of 5TeV/c in the central rapidity range and out to 1.5TeV/c in the 
forward rapidity region. Figure EHl shows the predictions with uncertainty ranges for three 
rapidity intervals. Figure |^ shows the ratios of the 40 eigenvector basis set predictions 
to the central prediction. The cross section uncertainties near the kinematic limit at the 
LHC due to PDF uncertainties are similar in magnitude to those obtained for Run 2 at the 
Tevatron. Again the extremes of the predictions are provided by the eigenvector basis sets 
+15 and —15, which correspond to extremes of the gluon distribution at large x. 

The jet cross section at the LHC will be a discovery mode for new physics. Supersymmetry 
would produce an enhanced cross section from production of gluino and squark jets. Extra 
dimensions would show up from production of Kaluza-Klein modes of the graviton. So 
the pure QCD prediction will be compared to the data and any difference will be the first 
evidence for these or other new physics scenarios. 
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FIG. 30: The inclusive jet cross section as a function of pT for three rapidity bins at the LHC. 
The three rapidity ranges are (0,1), (1,2) and (2,3). Predictions of all 40 eigenvector basis sets 
are superimposed. 
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FIG. 31: The uncertainty range of the inclusive jet cross section at the LHC. The curves are 
graphs of the ratios of the cross sections for the 40 eigenvector basis sets compared to the central 
(CTEQ6.1M) prediction (ordinate) versus p-r in GeV (ordinate). 
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VII. CONCLUSIONS 



Jet production serves as a probe of the highest scales accessible at a hadron-hadron 
collider. NLO QCD provides a good description of the inclusive jet cross sections measured 
by CDF and DO in Run 1 at the Tevatron, provided that the gluon distribution is enhanced 
at high X relative to earlier estimates. Such an outcome is now the standard consequence 
of the global fitting procedure due to the statistical power of the Tevatron jet data, and the 
preference for a larger gluon in all rapidity regions. The uncertainties on the jet cross section 
predictions still remain large, however, primarily due to the remaining uncertainties on the 
gluon distribution. In most regions, they are comparable to the experimental systematic 
uncertainties. 

The standard choice for the renormalization and factorization scales in evaluating inclu- 
sive jet cross sections has been E?^^ /2. We have re-done the global fitting using alternate 
scales of E^r^^ and 2E^r^^ and have found changes to the gluon distribution that are well within 
the gluon PDF uncertainty band. The effects of the resummation of threshold logarithms 
has been calculated for the first time in the forward rapidity regions. As for the central 
rapidity region, the effects of the threshold logarithms are relatively small. 

The increase of the center-of-mass energy to 1.96 TeV (from 1.8 TeV) and the expected 
greater integrated luminosity will lead to an expanded kinematic range for Run 2 at the 
Tevatron. Predictions have been made for both the inclusive jet cross section and its uncer- 
tainty for Run 2 and for the LHC. The PDF uncertainty on the jet cross section near the 
limits of the data reach in both cases are similar to those obtained for Run 1. 

We have carried out an exercise where we have examined the allowed range for any new 
physics such as compositeness, given the current scale of uncertainties on the jet cross section 
predictions. Given the Run lb jet data, any new physics must have a scale larger than 1.6 
TeV. When data is available from Run 2 at the Tevatron, and from the LHC, analogous 
analyses will be used to identify evidence for new physics or place even stricter limits on 
new physics. 
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Appendix 

One component of the theoretical uncertainty of any perturbative QCD calculation is the 
dependence of the answer on the various renormalization and factorization scales. One virtue 
of using NLO calculations in global determinations of parton distributions is the reduction 
in scale dependence usually exhibited by such calculations. The following brief discussion 
outlines the origin of the reduction of the scale dependence which is often observed with 
NLO calculations. 

Consider a large transverse momentum process such as the single jet inclusive cross section 
involving only massless partons. Furthermore, in order to simplify the notation, suppose 
that the transverse momentum is sufficiently large that only the quark distributions need 
be considered. In the following, a sum over quark flavors is implied. Schematically, one can 
write the lowest order cross section as 

E— = o = a^{^)aB®q{M)®q{M) (1) 

where a(/i) = as{^Ji)/'2'K and the lowest order parton-parton scattering cross section is de- 
noted hy as- The renormalization and factorization scales are denoted by /i and M, respec- 
tively. In addition, various overall factors have been absorbed into the definition of 0"^. The 
symbol ® denotes a convolution defined as 

f^9= f-f{-)9{y)- (2) 
Jx y y 

When one calculates the 0{a^) contributions to the inclusive cross section, the result can 
be written as 

a = a'^{^)aB®q{M)®q{M) 

+ 2a^(/i)61n(/i/pT)'7ij®g(M) (g)g(M) 

+ 2a3(^) ln(pr/M)Pg, ® ® g(M) ® q{M) 

+ a^{iJ,)K ®q{M)®q{M). (3) 

In writing Eq. Q, specific logarithms associated with the running coupling and the scale 
dependence of the parton distributions have been explicitly displayed; the remaining higher 
order corrections have been collected in the function K in the last line of Eq. (jHl). The n 
dependence of the running coupling is obtained from 

= /?(«(/^)) (4) 
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where (3 = — 6a^(l + ca) with 

, 33 - 2/ 153 - 19/ 

= and c = 5 

6 2(33 -2/) ^ ^ 

where / denotes the number of flavors. The scale dependence of the quark distributions is 
given by the nonsinglet DGLAP equation 

M^^|j^ = a(M)P,,®g(M). (6) 

Now, the NLO expression in Eq. Q should be independent of the choice of scale up to 
corrections of order a^. That this is indeed the case is easy to see. For convenience, set 
M = fi and calculate the logarithmic derivative of cr with respect to /i, 

^iTr- = 2a(/i)[-6a^(/i)(l + ca{^))]a ® q{fi) (g) 
a/i 

-2a^{n)Pgg (g) o-fi ® H 

= + O(a4). (7) 

The derivative of the lowest order contribution yields the first two lines of Eq. ((7j) while 
the derivative of the logarithms in the second and third lines of Eq. Q give rise to the last 
two lines in Eq. ((Tj) which cancel the fi dependent terms coming from the derivative of the 
lowest-order contribution. All other derivatives yield contributions which are of order a'^. 
Thus, the structure of the NLO calculation is such that the scale dependence is cancelled 
up to terms of order a^. 

The preceding discussion shows that one might expect a reduced scale dependence for 
the NLO calculation as opposed to that observed at lowest order. However, the uncancelled 
0{a^) and higher terms will still result in some scale dependence. It is easy to understand 
the systematics of the remaining dependence for the case at hand. Consider, first, the lowest- 
order contribution shown in the first line of Eq. (jS)). At large values of pt as xt = 2pj'/y/s 
tends to one, the scaling violations of the parton distributions lead to a decrease in the 
cross section as the scale increases. In addition, the running coupling decreases as the scale 
is increased. Therefore, at fixed px the lowest order contribution exhibits a monotonically 
decreasing behavior as the scale increases. 
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The second line of Eq. Q contains the contribution which partially offsets the scale 
dependence of the running coupling. Apart from the explicit logarithm, this contribution 
has much the same scale dependence as the lowest order contribution. Thus, for fj, < px this 
term is negative, becoming positive as the scale is increased beyond pt- 

The third line contains the contribution which partially offsets the scale dependence 
originating in the parton distributions. In the kinematic region under consideration, the 
convolution of the splitting function Pqq with the lowest contribution gives a negative result, 
reflecting the fact that the parton distributions decrease with increasing scale at large values 
of X. Therefore, this contribution is negative when the factorization scale M is chosen to be 
less than p^. 

Finally, the remaining term on line four of Eq. (j2)) has much the same scale dependence 
as the lowest order term, as the dominant scale dependence has already been exhibited in 
the terms on lines two and three. Thus, lines one and four give contributions which decrease 
monotonically with increasing scale while lines two and three start out negative, reach zero 
when the scales are equal to pt, and are positive for larger scales. This behavior explains why 
the typical NLO scale dependence in the high-p^ region yields a result which first increases 
with increasing scale, reaches a maximum, and then decreases. For many high-p^" processes 
involving massless partons this maximum occurs near fi = M = px- 

Next, it is interesting to consider how the location of the maximum changes with p^- 
Suppose that one goes to small values of xt such that the typical values of x in the parton 
distributions is around .2 or smaller. In this region the scaling violations change sign, 
corresponding to a growth of the parton distributions with increasing scale, rather than a 
decrease. The contribution of the third line of Eq.(l) then starts to compete against that 
of the second line. Accordingly, one must go to smaller values of /i/pr before seeing the 
turnover and decrease with decreasing scale. That is, the position of the maximum at fixed 
Pt moves to lower values of /i/pr as px decreases. The converse is also clearly true; as pt 
increases the position of the maximum moves to higher values of fi/pr- Thus, if one were 
to try to find an optimized scale at each value of pr, the pattern that would emerge would 
be that the optimal scale would be a fractional multiple of pt at lower px values and would 
increase to a multiple of pt greater than one at higher pt values. This pattern is precisely 
what is observed when using the single scale version of the Principle of Minimal Sensitivity 

Q 
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Another manifestation of this pattern concerns the behavior of the "K factor" , defined 
as the ratio of the NLO to LO result. At a fixed value of pr, the K factor will first rise and 
then fall with increasing scale. Suppose that at some pt value the maximum occurs for the 
choice n = M = Pt /2. Then, as one goes to higher p^ values the previous discussion shows 
that the maximum will move to a higher scale choice such as pt or even 2pT- At the larger 
Pt values it will then appear that the K factor may increase with increasing scale. This is 
just the behavior seen in the calculations done for the high rapidity DO jet data. 
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